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ABSTRACT 



Context. The possibility of observationally discriminating between various types of neutron stars, described by different equations of 
state of the nuclear matter, as well as differentiating neutron stars from other types of exotic objects, like, for example, quark stars, is 
one of the fundamental problems in contemporary astrophysics. 

Aims. We consider and investigate carefully the possibility that different types of rapidly rotating neutron stars, as well as other type 
of compact general relativistic objects, can be differentiated from the study of the emission properties of the accretion disks around 
them. 

Methods. We obtain the energy flux, the temperature distribution and the emission spectrum from the accretion disks around several 
classes of rapidly rotating neutron stars, described by different equations of state of the neutron matter, and for quark stars, described 
by the MIT bag model equation of state and in the CFL (Color-Flavor-Locked) phase, respectively. 

Results. Particular signatures appear in the electromagnetic spectrum, thus leading to the possibility of directly testing the equation 
of state of the dense matter by using astrophysical observations of the emission spectra from accretion disks. 

Key words, neutron stars - quark stars - accretion disks 



1. Introduction 

The quark structure of the nucleon indicates the possibility of a 
phase transition of confined hadronic matter to absolutely stable 
' strange quark matter at high densities (Itoh 1970 Bodmer 1971 
I IWitten 19841 1, which is referred to as strange matter hypothe- 
sis. If the hypothesis is true, then some neutron stars should 
actually be stars made of strange quark matter (strange stars) 
dAlcock et al. 19861 IHaensel et al. 1986l >. For a general review 
of strange star properties see Cheng et al. (1998a). 

There are several proposed mechanisms for the formation 
of quark stars. Quark stars are expected to form during the col- 
lapse of the core of a massive star, after the supernova explo- 
sion, as a result of a first or second order phase transition, re- 
sulting in deconfined quark matter dDai et al. 19951 1. The proto- 
neutron star core or the neutron star core is a favorable environ- 
ment for the conversion of ordinary matter to strange quark mat- 
ter dCheng et al. 1998b|IChan et al. 20091 . Another possibility is 
that some neutron stars in low-mass X-ray binaries can accrete 
sufficient mass to undergo a phase transition to become strange 
stars (Cheng & Dai 1996] l. This mechanism has also been pro- 
posed as a source of radiation emission for cosmological y-ray 
bursts ( |Cheng & Dai 1998) . 

Based on numerical integration of the general relativistic 
hydrostatic equilibrium equations a complete description of 
the basic astrophysical properties (mass, radius, eccentricity, 
Keplerian frequency etc.) of both static and rotating strange stars 
has been obtained, for different values of the bag constant and 
for different equations of state of the strange star (IWitten 19841 

IHaensel et al. 19861 IGondek-Rosinska et al. 2"000l 

|Dey et al. 1998] |Harko & Cheng 2002| i. Rotational proper- 



ties can discriminate between neutron and quark stars. Strange 
stars can reach much shorter periods than neutron stars, of 
the order of 0.5 ms ( |Cheng et al. 1998a i. r-mode instabilities 
in rapidly rotating strange stars lead to specific signatures in 
the evolution of pulsars with periods below 2.5 ms, and some 
data on pulsar properties are consistent with this assumption 
(Madsen 2000). Strange stars could have a radius significantly 
less than that of neutron stars (Cheng et al. 1998a I. 



Photon emissivity is the basic parameter for determin- 
ing macroscopic properties of stellar type objects. Because 
of the very high plasma frequency co p near the strange mat- 
ter edge, photon emissivity of strange matter is very low 
dAlcock et al. 1986b . The spectrum of equilibrium photons is 
very hard, with fra> > 20 MeV. The problem of the soft photon 
emissivity of quark matter at the surface of strange stars has also 
been considered dCheng & Harko 2003| [Harko & Cheng 2005 1. 
By taking into account the Landau-Pomeranchuk-Migdal effect 
and the absorption of the radiation in the external electron layer, 
the emissivity of the quark matter can be six orders of magnitude 
lower than the equilibrium black body radiation. 



The Coulomb barrier at the quark surface of a hot strange star 
may also be a powerful source of <? + e~ pairs, which are created 
in the extremely strong electric field of the barrier. At surface 
temperatures of around 10 11 K, the luminosity of the outflow- 
ing plasma may be of the order ~ 10 51 ergs -1 ([Usov 1998a 
lUsov 1998bl |Harko & Cheng 2006[ ). Moreover, for about one 
day for normal quark matter and for up to a hundred years for 
superconducting quark matter, the thermal luminosity from the 
star's surface, due to both photon emission and e + e~ pair pro- 
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duction, may be orders of magnitude higher than the Eddington 
limit dPage & Usov 2002) . 

However, despite these very specific signatures for quark 
stars, a definite method for discriminating them with respect to 
the neutron stars is still missing. 

It is generally expected that most of the astrophysical objects 
grow substantially in mass via accretion. Accretion discs are flat- 
tened astronomical objects made of rapidly rotating gas which 
slowly spirals onto a central gravitating body, with its gravita- 
tional energy degraded to heat. A fraction of the heat converts 
into radiation, which partially escapes, and cools down the ac- 
cretion disc. The only information that we have about accretion 
disk physics comes from this radiation, when it reaches radio, 
optical and X-ray telescopes, allowing astronomers to analyze 
its electromagnetic spectrum, and its time variability. The ef- 
ficient cooling via the radiation over the disk surface prevents 
the disk from cumulating the heat generated by stresses and dy- 
namical friction. In turn, this equilibrium causes the disk to sta- 
bilize its thin vertical size. The thin disk has an inner edge at 
the marginally stable orbit of the compact object potential, and 
the accreting plasma has a Keplerian motion in higher orbits 
dPage & Thorne 1974||Thbrne 1974). 



Recent astrophysical observations have shown that around 
many compact objects, as well as around black hole candidates, 
there are gas clouds surrounding the central compact object, and 
an associated accretion disc, on a variety of scales that could 
range up to a tenth of a parsec or even to a few hundred parsecs 
( Urry & Padovani 1995 ). These clouds are assumed to form a 
geometrically and optically thick torus (or warped disc), which 
absorbs most of the ultraviolet radiation and the soft X-rays. The 
gas exists in either the molecular or the atomic phase. Hence, 
important astrophysical information can be obtained from the 
observation of the motion of the gas streams in the gravitational 
field of compact objects. 

The first comprehensive theory of accretion disks was con- 
structed in Shakura & Sunyaev (1973) by using a Newtonian 
approach. This theory was extended to the general relativis- 
tic models of the mass accretion onto rotating black holes 
dNovikov & Thorne 1973K These pioneering works developed 
thin steady-state accretion disks, where the accreting matter 
moves in Keplerian orbits. The hydrodynamical equilibrium in 
the disk is maintained by an efficient cooling mechanism via ra- 
diation transport. The photon flux emitted by the disk surface 
was studied under the assumption that the disk emits a black 
body radiation. The properties of radiant energy flux over the 
thin accretion disks were further analyzed, with the effects of 
the photon capture by the hole on the spin evolution taken into 
account as well ([Page & Thorne 1974 IThorne 1974K In these 



works the efficiency with which black holes convert rest mass 
into outgoing radiation in the accretion process was also com- 
puted. 

The disc and boundary layer luminosity for accret- 
ing rapidly rotating neutron stars with low magnetic fields 
in a fully general relativistic manner were considered in 
Tha mpan & Datta (1998)| The effect of the quadrupole compo- 
nent in the mass distribution of a rapidly rotating neutron star 
on the energy release in the equatorial (or boundary) layer on 
the surface of the accreting star and in the accretion disk in 
the cases where the stellar radius is smaller (or larger) than 
the radius of the last stable circular orbit was analyzed in 
Sibgatullin & Sunyaev (1998)| The temperature profiles of (thin) 
accretion disks around rapidly rotating neutron stars (with low 
surface magnetic fields), taking into account the full effects of 
general relativity were obtained in Bhattacharyya et al. (2000) 



The calculations suggest that the neutron star in Cygnus X-2 
rotates close to the centrifugal mass-shed limit. General rel- 
ativistic spectra from accretion disks around rotating neutron 
stars in the appropriate spacetime geometry for several differ- 
ent equations of state, spin rates, and masses of the compact ob- 
ject have been also computed (Bhattacharyya et al. 2001b). The 
spectra from the accretion disks around a rapidly rotating neu- 
tron star, by taking into account the Doppler shift, the gravi- 
tational redshift and the light-bending effects were obtained in 
Bhattacharyya et al. (2001) Light bending significantly modi- 



fies the high-energy part of the spectrum. These results could 
be important for modeling the observed X-ray spectra of low- 
mass X-ray binaries, containing fast-spinning neutron stars. The 
temperature profiles of accretion discs around rapidly rotating 
strange stars, using constant gravitational mass equilibrium se- 
quences of these objects, and considering the full effect of gen- 
eral relativity were computed in [Bha ttacharyy a et al. (2001a)] 
The structure parameters and temperature profiles obtained were 
compared with those of neutron stars, as an attempt to provide 
signatures for distinguishing between the two. The obtained re- 
sults implies the possibility of distinguishing these objects from 
each other by sensitive observations through future X-ray de- 
tectors. The accretion disc temperature profiles, disc luminosity 
and boundary layer luminosity for rapidly rotating neutron stars, 
considering the full effect of general relativity, were calculated 



in|Bhattacharyya (2002) 



The emissivity properties of the accretion disks have been 
recently investigated for a large class of compact objects, 
such as rotating and non-rotating boson or fermion stars 
dTorres 20021 1 Yuan et al. 20041 IGuzman 2006b, as well as for 
the modified f(R) type theories of gravity dPun et al. 2008al l, for 
brane world black holes dPun et al . 2008b ), and for wormholes 
dHarko et al. 20 08 Harko et al 2009). The radiation power per 
unit area, the temperature of the disk and the spectrum of the 
emitted radiation were given, and compared with the case of a 
Schwarzschild or Kerr-Newman black holes of an equal mass. 

It is the purpose of the present paper to consider a compara- 
tive systematic study of the properties of the thin accretion disks 
around rapidly rotating neutron and strange stars, respectively, 
and to obtain the basic physical parameters describing the disk, 
like the emitted energy flux, the temperature distribution on the 
surface of the disk, as well as the spectrum of the emitted equi- 
librium black body radiation. 

In order to obtain the emissivity properties of the disk, the 
metric outside the rotating general relativistic stars must be de- 
termined. In the present study we calculate the equilibrium con- 
figurations of the rotating neutron and quark stars by using the 
RNS code, as introduc ed in[St ergioulas & Friedman (1995)] and 
discussed in detail in Stergioulas et al. (2003)| This code was 
used for the study of different models of rotating neutron stars 
in Nozawaet al. (1998) and for the study of the rapidly rotat- 
ing strange stars ( |Stergioulas et al. 1999| l. The software provides 
the metric potentials for various types of compact rotating gen- 
eral relativistic objects, which can be used to obtain the physical 
properties of the accretion disks. Particular signatures appear in 
the electromagnetic spectrum, thus leading to the possibility of 
directly testing the equation of state of the dense matter by using 
astrophysical observations of the emission spectra from accre- 
tion disks. 

The present paper is organized as follows. The properties of 
the general relativistic thin accretion disks onto compact objects 
are briefly described in Section II. The equations of state used 
in the present study are presented in Section III. In Section IV 
we consider the radiation flux, spectrum and efficiency of thin 
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accretion disks onto several classes of neutron stars and quark 
stars. We discuss and conclude our results in Section V. 



2. Thin accretion disks onto general relativistic 
compact objects 

A thin accretion disk is an accretion disk such that in cylindri- 
cal coordinates (r, <p, z) most of the matter lies close to the ra- 
dial plane. For the thin accretion disk its vertical size (defined 
along the z-axis) is negligible, as compared to its horizontal ex- 
tension (defined along the radial direction r), i.e, the disk height 
H, equal to the maximum half thickness of the disk, is always 
much smaller than the characteristic radius R of the disk, H <K R 
(Novi kov & Thorne 1973I I. The thin disk is in hydrodynamical 
equilibrium, and the pressure gradient and a vertical entropy gra- 
dient in the accreting matter are negligible. In the steady-state 
accretion disk models, the mass accretion rate Mo is supposed 
to be constant in time, and the physical quantities of the accret- 
ing matter are averaged over a characteristic time scale, e.g. At, 
and over the azimuthal angle A<p = 2n, for a total period of the 
orbits and for the height H. The plasma moves in Keplerian or- 
bits around the compact object, with a rotational velocity Q.. The 
plasma particles have a specific energy E, and a specific angu- 
lar momentum L, which depend only on the radii of the orbits. 
The four-velocity of the particles is vP, and the disk has an av- 
eraged surface density Z. The accreted matter is modeled by an 
anisotropic fluid source, where the density po (the specific heat 
is neglected), the energy flow vector <f, and the stress tensor f* 
are measured in the averaged rest-frame. The energy-momentum 
tensor takes the form (Novikov & Thor ne 19731 1 



(1) 



where = 0, M /J f /iv = 0. The four-vectors of the energy and 
of the angular momentum flux are defined by —W = P , v (<9/<9?) v 
and J M = T^yid/dcpY, respectively. The four dimensional con- 
servation laws of the rest mass, of the energy and of the angular 
momentum of the plasma provide the structure equations of the 
thin disk. By integrating the equation of the rest mass conserva- 
tion, Vfiipou 11 ) = 0, it follows that the time averaged accretion 
rate Mq is independent of the disk radius: 



Mo = -2nrLu r = const , 



where a dot represents the derivative with respect to the time co- 
ordinate (Page & Thorne 1974b. The averaged rest mass density 

TR 

is defined by S(r) = J _ H (po)dz, where (po) is the rest mass den- 
sity averaged over At and 2n. The conservation law V ll E ,i — of 
the energy can be written in an integral form as 



[M E - InrQW/Xr = 4tt yf^gFE 



(3) 



where a comma denotes the derivative with respect to the ra- 
dial coordinate r. Eq. (01 shows the balance between the energy 
transported by the rest mass flow, the dynamical stresses in the 
disk, and the energy radiated away from the surface of the disk, 
respectively. The torque W/ in Eq. (O is given by 



, r = f (t/)dz, 

J-H 



where (f/ ) is the <j> — r component of the stress tensor, averaged 
over At and over a 2n angle. The law of the angular momentum 
conservation, V^/ 1 " = 0, states in its integral form the balance of 
the three forms of the angular momentum transport, 



[M Q L - 2nrW/], r = 4tt iJ=gFL 



By eliminating W/ from Eqs. Q and (0, and by applying 
the universal energy -angular momentum relation dE = QdJ for 
circular geodesic orbits in the form E t r - ClL r , the flux of the 
radiant energy over the disk can be expressed in terms of the 
specific energy, angular momentum and the angular velocity of 
the black hole. Then the flux integral leads to the expression of 
the energy flux F(r), which is given by ( |Page & Thorne 1974| 
IThorne 19741 



F(r) 



Mo 



4tf V=? (E - QL) 



E-QL) 2 J,„„ 



(E - QL)L ?r dr . 



(6) 



where the no-torque inner boundary conditions were also pre- 
scribed (Page & Thorne 1974b . This means that the torque van- 
ishes at the inner edge of the disk, since the matter at the 
marginally stable orbit r ms falls freely into the black hole, and 
cannot exert considerable torque on the disk. The latter assump- 
tion is valid as long as strong magnetic fields do not exist in the 
plunging region, where matter falls into the hole. 

Once the geometry of the space-time is known, we can de- 
rive the time averaged radial distribution of photon emission for 
accretion disks around black holes, and determine the efficiency 
of conversion of the rest mass into outgoing radiation. After ob- 
taining the radial dependence of the angular velocity Q., of the 
specific energy E and of the specific angular momentum L of 
the particles moving on circular orbits around the black holes, 
respectively, we can compute the flux integral given by Eq. ©. 

Let us consider an arbitrary stationary and axially symmetric 
geometry, 



ds — g, t dt + : 



+ g rr dr 2 + g ge d6 2 + g^dcf) 2 



(1) 



where in the equatorial approximation (\9 - n/2\ <K l)the met- 
ric functions g tt , g,^, g rr , ggg and g^ depend only on the radial 
coordinate r. The geodesic equations take the form 



dt _ Egtf + Lg t4 , d<p _ Eg, $ + Lg„ 
dr g 1 - s„Sd,d, ' dr 

and 



,2) r \'£\ =V(r), 



(8) 



(9) 



respectively, where r is the affine parameter, and the potential 
term V(r) is defined by 



V(r) = -l + 



F 2 g^ + 2ELg t4 , + L 2 g„ 



(10) 



For circular orbits in the equatorial plane the conditions 
V(r) = and V r (r) = 0, respectively, must hold for all disk 
configurations. From these conditions one can obtain the spe- 
cific energy E, the specific angular momentum L and the angu- 
lar velocity O of the particles moving on circular orbits around 
spinning general relativistic stars as 



(4) =r 



E - — 



gtt + gt^ 



L = 



y/-gn ~ 2gt^ ~ g<p<pQ 2 
yj-gtt ~ 2g l<t> n - g^Q. 2 



(5) 



„ _ d4 _ -gt<j,,r+ V(gf«A.r) 2 ~ gtt.rg 
dt gtbtb.r 



(ID 

(12) 
(13) 
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The marginally stable orbits r ms around the central object are 
determined by the condition V„ .(r)| ( _ = 0, which is equivalent 
to the equation 

[E 2 g^ rr + 2ELg,^ rr + L 2 g, Urr - (gy - 

By inserting Eqs. (fTTTi-(fT2li into Eq. ( fT4l i. and solving the 
resulting equation for r ms , we obtain the marginally stable orbits, 
once the metric coefficients g„, g t ^ and g^ are explicitly given. 

The accreting matter in the steady-state thin disk model is 
supposed to be in thermodynamical equilibrium. Therefore the 
radiation emitted by the disk surface can be considered as a per- 
fect black body radiation, where the energy flux is given by 
F(r) = crT 4 (r) (cr is the Stefan-Boltzmann constant), and the 
observed luminosity L (v) has a redshifted black body spectrum 
dTorres 2002} : 



L (v) = And I (v) = — r- cos y 

TIC 



r r f r 2rr 

Jr, JO ' 



v 3 e rd(f>dr 



exp(hvJT) - 1 



(15) 



Here d is the distance to the source, I(v) is the Planck dis- 
tribution function, y is the disk inclination angle, and r, and rj 
indicate the position of the inner and outer edge of the disk, re- 
spectively. We take r, = r ms and rj — » oo, since we expect the 
flux over the disk surface vanishes at r — » oo for any kind of gen- 
eral relativistic compact object geometry. The emitted frequency 
is given by v e — v(l + z), where the redshift factor can be written 
as 



l+Z 



1 + Or sin sin y 
y/-g,t - 2Qg ¥ - Q. 2 gt 



(16) 



where we have neglected the light bending (ILuminet 19791 
|Bhattacharyya et al. 2001b) . 

The flux and the emission spectrum of the accretion disks 
around compact objects satisfy some simple scaling relations, 
with respect to the simple scaling transformation of the radial 
coordinate, given by r — > T = r/M, where M is the mass of the 
black hole. Generally, the metric tensor coefficients are invariant 
with respect of this transformation, while the specific energy, the 
angular momentum and the angular velocity transform as E — > 
E, L — > ML and Q. — > Q./M, respectively. The flux scales as 
F(r) — > F(F)/M 4 , giving the simple transformation law of the 
temperature as T(r) — > T (r) /M. By also rescaling the frequency 
of the emitted radiation as v — > ~v — v/M, the luminosity of the 
disk is given by L (v) — » L (v) /M. On the other hand, the flux is 
proportional to the accretion rate Mo, and therefore an increase 
in the accretion rate leads to a linear increase in the radiation 
emission flux from the disc. 

The efficiency e with which the central object converts rest 
mass into outgoing radiation is the other important physical pa- 
rameter characterizing the properties of the accretion disks. The 
efficiency is defined by the ratio of two rates measured at in- 
finity: the rate of the radiation of the energy of the photons es- 
caping from the disk surface to infinity, and the rate at which 
mass-energy is transported to the compact object. If all the emit- 
ted photons can escape to infinity, the efficiency depends only on 
the specific energy measured at the marginally stable orbit r ms , 



1 



\r=r ms 



(17) 



There are some solutions for the geometry of the neutron and 
quark stars where r ms is at the surface of the star or even takes 
values less than the surface radius of the central object. In this 
case the inner edge of the disk touches the surface of the star 



and plasma under the effect of any perturbation due to hydro- 
or magnetohydrodynamic instabilities in the disk will leave the 
disk and hit the surface. Then the energy E e transferred to the 
star from the disk is measured at the radius R e of the star, and 
the efficiency takes the form 



e = 1 - E e 



(18) 



where E e = £j , For Schwarzschild black holes the efficiency 
is about 6%, no matter if we consider the photon capture by 
the black hole, or not. Ignoring the capture of radiation by the 
black hole, e is found to be 42% for rapidly rotating black holes, 
whereas the efficiency is 40% with photon capture in the Kerr 
potential. For neutron and quark stars the efficiency is varying in 
a broader range, especially if we take into account that E ms and 
E e can have very different values for different neutron and quark 
star models. 



3. Equations of state 

In order to obtain a consistent and realistic physical description 
of the rotating general relativistic neutron and quark stars, as a 
first step we have to adopt the equations of state for the dense 
neutron and quark matter, respectively. In the present study of 
the accretion disks onto rapidly rotating neutron and quark stars 
we consider the following equations of state: 

1) Akmal-P andharipande-Ravenhall (APR) EOS 
(lAkmaletal. 19981 1. EOS APR has been obtained by using 
the variational chain summation methods and the Argonne vis 
two-nucleon interaction. Boost corrections to the two-nucleon 
interaction, which give the leading relativistic effect of order 
(v/c) 2 , as well as three-nucleon interactions, are also included 
in the nuclear Hamiltonian. The density range is from 2 x 10 14 
g/cm 3 to 2.6 x 10 15 g/cm 3 . The maximum mass limit in the static 
case for this EOS is 2.2OM . We join this equation of state to the 



composite BBP (e/c 2 > 4.3 x 10 n g/cm 3 ) (IBaymetal. 1971a i 

, , n A 1 ._ -—11 1 . I ' - , I -L 



BPS (10 4 g/cm 3 < 4.3 x 10 n g/cm 3 ) (IBaymetal. 1971b) - FMT 



(e/c 2 < 10 4 g/cm 3 ) ( Feynman et al. 1949) equations of state 



respectively. 

2) Douchin-Haensel (DH) EOS dDouchin & Haensel 200Tb . 
EOS DH is an equation of state of the neutron star matter, de- 
scribing both the neutron star crust and the liquid core. It is 
based on the effective nuclear interaction SLy of the Skyrme 
type, which is particularly suitable for the application to the cal- 
culation of the properties of very neutron rich matter. The struc- 
ture of the crust, and its EOS, is calculated in the zero temper- 
ature approximation, and under the assumption of the ground 
state composition. The EOS of the liquid core is calculated as- 
suming (minimal) npe/j composition. The density range is from 
3.49 x 10" g/cm 3 to 4.04 x 10 15 g/cm 3 . The minimum and maxi- 
mum masses of the static neutron stars for this EOS are 0.094M o 
and 2.05 M , respectively. 

3) Shen-To ki-Oyamatsu-Sumiyoshi (STOS) EOS 
(She n et al. 1998K The STOS equation of state of nuclear 
matter is obtained by using the relativistic mean field theory 
with nonlinear cr and co terms in a wide density and temperature 
range, with various proton fractions. The EOS was specifically 
designed for the use of supernova simulation and for the neutron 
star calculations. The Thomas-Fermi approximation is used to 
describe inhomogeneous matter, where heavy nuclei are formed 
together with free nucleon gas. We consider the STOS EOS 
for several temperatures, namely T = 0, T — 0.5 and T = 1 
MeV, respectively. The temperature is mentioned for each STOS 
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equation of state, so that, for example, STOS represents the 
STOS EOS for T — 0. For the proton fraction we chose the 
value Y p = 0.001 in order to avoid the negative pressure regime 
for low baryon mass densities. 

4) Relativistic Mean Field (RMF) equations of state with 
isovector scalar mean field corresponding to the (5-meson- RMF 
soft and RMF stiff EOS (IKubis & Kutschera 1997) . While the 
(5-meson mean field vanishes in symmetric nuclear matter, it 
can influence properties of asymmetric nuclear matter in neu- 
tron stars. The Relativistic mean field contribution due to the 
(5-field to the nuclear symmetry energy is negative. The energy 
per particle of neutron matter is then larger at high densities than 
the one with no <5-field included. Also, the proton fraction of B- 
stable matter increases. Splitting of proton and neutron effective 
masses due to the (5-field can affect transport properties of neu- 
tron star matter. The equations of state can be parameterized by 
the coupling parameters = gl/m 2 , C 2 = g 2 Jml, b = bjg\ 
and c = c/g 4 , where m a and m w are the masses of the respective 
mesons, and b and c are the coefficients in the potential energy 
U (cr) of the <x-field. The soft RMF EOS is parameterized by 
C 2 = 1.582 fm 2 , C 2 = 1.019 fm 2 , b = -0.7188 and c = 6.563, 
while the stiff RMF EOS is parameterized by C 2 = 1 1.25 fm 2 , 
C 2 = 6.483 fm 2 , b = 0.003825 and c = 3.5 x 10' 6 , respectively . 

5) Baldo-Bombaci-Burgio (BBB) EOS dBaldo et al. 79971 ) . 
The BBB EOS is an EOS for asymmetric nuclear matter, de- 
rived from the Brueckner-Bethe-Goldstone many-body theory 
with explicit three-body forces. Two EOS's are obtained, one 
corresponding to the Argonne AV14 (BBBAV14), and the other 
to the Paris two-body nuclear force (BBBParis), implemented by 
the Urbana model for the three-body force. The maximum static 
mass configurations are M max - 1.8M and M max = 1.94M 
when the AVI 4 and Paris interactions are used, respectively. The 
onset of direct Urea processes occurs at densities n > 0.65 fm" 3 
for the AV14 potential and n > 0.54 firT 3 for the Paris po- 
tential. The comparison with other microscopic models for the 
EOS shows noticeable differences. The density range is from 
1.35 x 10 14 g/cm 3 to 3.507 x 10 15 g/cm 3 . 

6) Bag model equation of state (Q) EOS dWitten 19841 
Cheng et al. 1998a) . For the description of the quark matter we 
adopt first a simple phenomenological description, based on the 
MIT bag model equation of state, in which the pressure p is re- 
lated to the energy density p by 



10' 



,36 



p = i(p-4fi), 



(19) 



where B is the difference between the energy density of the per- 
turbative and non-perturbative QCD vacuum (the bag constant), 
with the value B = 4.2 x 10 14 g/cm 3 . 

7) Color-Flavor-Locked (CFL) EOS It is generally agreed 
today that the color-flavor-locked (CFL) state is likely to be 
the ground state of matter, at least for asymptotic densities, 
and even if the quark masses are unequal dAlford et al. 1999, 
Rapp et al. 2000| lAlford et al. 20071). Moreover, the equal num- 



ber of flavors is enforced by symmetry, and electrons are absent, 
since the mixture is automatically neutral. By assuming that the 
mass m s of the s quark is not large as compared to the chemical 
potential p, the thermodynamical potential of the quark matter in 
CFL phase can be approximated as (Lugones & Horvath 2002) 



CFL - 



3p 3m 2 

47T 2 47T 2 



1 - 121n(m s /2;u) 4 3 2 2 

_ ms --Ap + B,(20) 



where A is the gap energy. With the use of this expression the 
pressure P of the quark matter in the CFL phase can be ob- 
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Fig. 1. Pressure as a function of density (in a logarithmic scale) 
for the equations of state Q, CFL 150, CFL300, DH, RMF soft, 
RMF stiff, STOS 0. STOS 0.5, STOS 1, BBBAV14, BBBParis 
and APR, respectively. 



tained as an explicit function of the energy density s in the form 
( jLugones & Horvath 2002] ) 



P = - (e 
3 v 

where 



45)- 



2A 2 6 2 m 2 s d 2 
^ 2_ + ^ 2 "' 



,r =-a+ -\\a 2 + -n 2 (s-B), 



(21) 



(22) 



and a = -m 2 /6 + 2A 2 /3. In the following the value of 
the gap energy A considered in each case will be also 
mentioned for the CFL equation of state, so that, for ex- 
ample, CFL 150 represents the CFL EOS with A = 150 
MeV. For the mass of the strange quark we adopt the value 
m s = 150 MeV. The maximum mass of the stars in the 
CFL model is given by M max = 1.96M Q (1 + 6) / VB~6o, 
with 6 = 0.15(A/100MeV) 2 (60MeVfirT 3 /fi) 
(Horvath & Lugone s 2004[ ). Hence large maximum masses 
(M > 4M ) can be found for standard values of B and 
A > 250 MeV. For a discussion of the maximum mass of 
both rotating and non-rotating general relativistic objects see 
Stergioulas et al. (2003) For rotating stars the maximum mass 

can be of the order of M max w 6.1M (2 x 10 14 g crrT 3 /p m ) 1/2 , 
where p m is the matching density. 

The pressure-density relation is presented for the considered 
equations of state in Fig.[TJ 



4. Thin disk accretion onto neutron and quark stars 

In the present Section we consider the properties of the general 
relativistic thin accretion disks formed around neutron stars and 
quark stars. The metric outside the rotating compact general rel- 
ativistic stars can be described, in quasi-isotropic coordinates, 
as 



ds 2 = -e7 + Pdt 2 +e m (dr 2 + r 2 dd 2 ) 



)+e f -Pf 2 sin 2 6(d(f>- 



udtf , (23) 



where the metric potentials y, p, a and the angular velocity of 
the stellar fluid relative to the local inertial frame a> are all func- 
tions of the quasi-isotropic radial coordinate f and of the po- 
lar angle Q. The RNS code produces the metric functions in a 
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quasi-spheroidal coordinate system, as functions of the parame- 
ter s - ?l (r + r e ), where r e is the equatorial radius of the star, 
which we have converted into Schwarzschild-type coordinates r 
according to the equation r = rexp [(y - p) /2]. To obtain the 
radius of the marginally (or innermost) stable circular orbits r ms 
we have used a truncated form of the analytical approximation 
given as dShibata & Sasaki 1998b , 

^ * (l - 0.54433^ - 0.2261V + 0.17989g 2 - 0.23002^+ 
0.26296^g 2 - 0.29693^ + 0.44546^ 2 2 2 ) , (24) 

where q = J/M 2 and Q2 = -M2/M 3 , respectively, and where 
J is the spin angular momentum, and Mi is the quadrupole mo- 
ment. Once the metric outside the rotating stars is known, from 
Eqs. (fT2b-([T3b we obtain the angular velocity, the specific energy 
and the specific angular momentum of the orbiting plasma parti- 
cles in the disk. Then we calculate the integral given by Eq. ©, 
measuring the photon flux emitted by the disk surface in ther- 
modynamical equilibrium. From the flux we obtain the temper- 
ature distribution of the disk, as well as the spectra of the emit- 
ted radiation. For all our calculations we use an accretion rate 
of M — 1 x 10~ 12 M o /yr, and assume that cosy = 1. In order 
to compare the thermal and emission properties of the thin ac- 
cretion disks onto neutron and quark stars, we consider several 
classes of models in which both the neutron and the quark stars 
have some common characteristics. 

4.1. Accretion disks onto neutron and quark stars with same 
mass and angular velocity 

As a first example of the electromagnetic radiation characteris- 
tics of the accretion disks we consider the properties of disks 
around rapidly rotating neutron and quark stars having the same 
gravitational mass M and angular velocity O. If M and Q. as well 
as the electromagnetic characteristics of the radiation emission 
from the disk could be determined by using astronomical and as- 
trophysical observations, one may attempt to identify the type of 
the central star by analyzing the differences in the physical pa- 
rameters (flux, temperature and spectrum) of the disk. The phys- 
ical parameters of the accretion disks are obtained for a fixed 
total mass and rotational frequency for all the equations of state 
discussed in Section[3] 

In Table Q] we present the main physical parameters of the 
considered neutron and quark stars, with the total mass set to 
values around 1.8M , and the angular velocity fixed to values 
around 5 x 10 3 s _1 . In this, and the following Tables, p c is the 
central density, M is the gravitational mass, Mo is the rest mass, 
R e is the circumferential radius at the equator, Q. is the angu- 
lar velocity, Q p is the angular velocity of a particle in circu- 
lar orbit at the equator, T/W is the rotational-gravitational en- 
ergy ratio, cJ/GM 2 is the angular momentum, / is the moment 
of inertia, O2 gives the mass quadrupole moment M2 so that 
M 2 = c 4 O 2 /(G 2 M 3 M ), h + is the height fr om the surface of the 
last stable co-rotating circular orbit in the equatorial plane, h- is 
the height from surface of the last stable counter-rotating circular 
orbit in the equatorial plane, u c /Q. is the ratio of the central value 
of the potential to to Q, r e is the coordinate equatorial radius, and 
r p /r e is the axes ratio (polar to equatorial), respectively. 

In Fig. [2] we present the total flux F (r) radiated by the accre- 
tion disk around the stellar models DH, RMF soft/stiff, STOS 
T = 0.0, 0.5 and 1.0, BBBAV14, BBBParis, APR, Q and 
CFL150, respectively, all the stars having the mass of 1.8M and 




Fig. 2. The time-averaged flux radiated by a thin accretion disk 
around rotating neutron and quark stars, with the same total mass 
M = 1.8M Q , and rotational velocity Q. - 5 x 10 3 s _1 . 



the angular velocity 5 x 10 3 s 1 . We can distinguish two types of 
behavior for the curves in this figure. 

The first group contains the models with relatively small stel- 
lar radii, of the order of 10 - 12 km. The stars with EOS APR, 
BBBAV14, BBBParis and DH belong to this group, and the in- 
ner edges of the accretion disk (indicated by the left boundary 
of the radial flux profile) are located at low dimensionless radii, 
close to each other, with r/M ^ 5.2-5.6. The shape of their flux 
profiles is rather similar, and only the stellar model APR pro- 
duces a maximal flux value higher by a factor of 1 .4 with respect 
to the other ones. The maximal flux values for the stars with 
the DH and BBBParis type EOS's are close to each other, and 
the model BBBAV14 has the smallest maximal value. For the 
models BBBAV14, BBBParis, and DH, respectively, the disk ra- 
diates its maximal flux at higher radii as compared to the APR 
EOS. The model RMF soft and the quark star model Q still be- 
long to this first group, and they have accretion disks with the 
inner edge at r/M ^ 5.4. Nevertheless, the flux profiles of their 
disk radiation are different. The highest values for F (r) among 
all the models of the first group are obtained for the RMF soft 
model. The flux emitted by the accretion disk around the quark 
star described by the bag model equation of state Q is lower by 
a factor of around 0.7, as compared to the fluxes emitted by the 
disks around the neutron stars, but the flux profile is very similar 
to the DH model profile. The only difference between them is in 
the location of the inner edge of the disk. The quark star model 
CFL150 also belongs to the first group. The inner disk edge is 
located at r/M ^5.6 for this type of EOS, but the maximal flux 
emitted by the disk is close to the one emitted by the BBBParis 
model. 

The RMF stiff and the STOS type neutron stars form a sec- 
ond group, since the accretion disks around them emit a flux that 
is much smaller, as compared to the fluxes of the first group. The 
stellar radii of these stars are much greater than those of the stars 
in the first group, and their marginally stable orbits are located 
at even higher radii, between r/M ^ 7 and 1 1 .7, respectively. In 
the STOS type models, the maximum values of their fluxes are 
increasing, and the radii of the left boundary of F(r) are decreas- 
ing, with increasing star temperatures. The EOS RMF stiff, with 
the smallest radius of the inner edge, gives the highest maxi- 
mal flux in the second group. The temperature distribution in 
the disk, shown in Fig. [3] follows the same pattern as the one 
described for the flux profiles. 
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Table 1. Physical parameters of the compact stars with total mass M as 1.8M and rotational frequency Q « 5 x 10 3 s 
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Fig. 3. The temperature distribution in a thin accretion disk 
around rotating neutron and quark stars with the same total mass 
M = I.SMq, and rotational velocity O = 5 x 10 3 s _1 . 
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Fig. 4. The emission spectrum of the accretion disk around ro- 
tating neutron and quark stars with the same total mass M = 
1.8M Q , and rotational velocity Q. - 5 x 10 3 s _1 . 



In Fig. [4] we present emission spectra of the disks, which 
similarly to the flux and temperature curves, also exhibit some 
distinctive features. Both the amplitude and the cut-off frequency 
of the spectra are determined roughly by the maximum values of 
F(r): for higher maximal flux values of an accretion disk around 
a given stellar model we obtain higher maximal amplitudes of 
the disk spectra, and they are shifted to higher frequencies. This 
is clearly seen in Fig. [4] for the neutron stars with the RMF stiff 
and the STOS type EOS's. However, there are some exceptions. 
The spectral maximum is lower for the model RMF soft as com- 
pared to the APR type EOS, but its cut-off is still located at a 
higher frequency. The spectral maxima are almost the same for 
the BBBAV14 and BBBParis models, but the cut-off frequency 
is higher for the BBBAV14 type EOS. The spectral maxima for 
all the models in the plot are located between 2xl0 16 and4xl0 16 
Hz. The spectral features of the neutron star model DH and of 
both types of quark star models are essentially the same, and thus 
only the morphological properties of the accretion disk, i.e., the 
radial distribution of the flux or the disk temperature, can be used 
to discriminate these type of compact stellar objects from each 
other. 



4.2. Accretion disks onto neutron and quark stars rotating 
with Keplerian velocities 

As a second case of comparison between the neutron and quark 
star properties we consider the Keplerian models of the neutron 
and quark stars, where they rotate at the maximal frequency al- 
lowed by their physical properties. We fix the central density for 
all the stellar models to p c = 10 15 g/cm 3 , and we compute the 
physical parameters of all the stellar models for the case of the 
Keplerian rotation. The results are represented in Table|2] 

The flux curves for these configurations are shown in Fig. [5] 
Here each EOS type has very distinctive features. Only the STOS 
type stars, and the models BBBAV14 and BBBParis produce 
similar flux profiles. Nevertheless, there is a temperature depen- 
dence for the models STOS in the location of the inner disk edge 
(r/M = 4.1 - 4.2), and of the maximal value of the photon flux 
emerging from the disk: with increasing temperatures we obtain 
greater radii for the inner edge of the disk, and the maximal flux 
values are inversely proportional to the temperature of the star. 
For the models BBBAV14 and BBBParis, we see that the loca- 
tion of the inner edge of the accretion disks around these types of 
neutron stars are essentially the same (r/M 7.1), but the model 
BBBAV14 produces a higher maximal flux than the BBBPAris 
type EOS does. Although the thin disk rotating around the APR 
type star has an inner edge at only a little bit lower radius than 
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Table 2. Physical parameters of the neutron and quark stars rotating at Keplerian frequencies. All stars have the same central density 

Pc- 



the disks for the models BBBAV14 and BBBParis, respectively, 
its maximal radiated away flux is a factor of 1.4 smaller than 
those for the neutron stars with the latter types of EOS. 

The inner edge of the accretion disk around the quark star 
CFL300 is located at the lowest radius (r/M ^ 3.6), as com- 
pared with the other EOS models in this configuration, whereas 
the model RMF soft has the greatest radius for the inner edge 
(r/M ^ 9.5). The flux values for the RMF soft type star are the 
smallest ones as well. The maximal flux for the model DH is the 
greatest in this group, and it is at least a factor of 3.5 greater than 
the one for the model RMS stiff. The flux profiles of the Q type 
quark star and of the neutron star with the RMF stiff type EOS 
are very similar, but the inner edge of the disk rotating around 
the neutron star is located at a lower radius (r/M ^ 5.5 for the 
quark star, and r/M ^5.1 for the neutron star). For the radial 
distributions of the disk temperature we obtain the same char- 
acteristic behavior and groups of models, which we present in 
Fig. [6] 

The spectral features, as seen in Fig. [7] are strongly influ- 
enced by the flux and the temperature distributions of the disk. 
For the case of the Keplerian rotation the emission spectra of 
the disk are rather different from those of the previous config- 
uration, with same gravitational mass and angular velocity. The 
proportionality between the location of the inner disk edge and 
the maximal amplitudes of the spectra holds for this case as 
well. Therefore the greatest spectral maximum is produced by 
the quark star model CFL300, whereas the smallest one is ob- 
tained for the model RMF stiff. The STOS models have also 
high amplitudes in their spectra, and their maximal values, in 
contrast with the previous case, are inversely proportional to the 
temperature. 

There are no big differences between the cut-off frequencies 
in the spectra, which are around 2 x 10 16 - 3 x 10 16 Hz. The 
spectrum belonging to the RMF stiff type EOS has its maxi- 
mum at the lowest frequency, and the spectral maxima for the 
models DH, BBBAV14 and BBBParis are located at the highest 
frequencies. Since the spectra for the latter two models are very 
similar, we can use only the differences in the flux maxima to 
discriminate between them. This is also the case for the quark 
star model Q and the RMF stiff type neutron star as well, where 
the only discriminating factor is the difference in the locations of 
the inner disk edges. The models APR and DH have also similar 
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Fig. 5. The time-averaged flux radiated by a thin accretion disk 
around neutron and quark stars rotating at the Keplerian velocity. 
All the stars have the same central density p c = 10 g/cm 3 . 
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Fig. 6. The temperature distribution by a thin accretion disk 
around neutron and quark stars rotating at the Keplerian velocity. 
All the stars have the same central density p c = 10 g/cm 3 . 



spectra, but both their maxima and cut-off frequencies exhibit 
some differences. 
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Fig. 7. The emission spectrum of the accretion disk around neu- 
tron and quark stars rotating at the Keplerian velocity. All the 
stars have the same central density p c = 10 15 g/cm 3 . 
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Fig. 9. Temperature distribution in a thin accretion disk around 
rotating neutron and quark stars, with the same central density 
(p c = 10 15 g/cm 3 ) and eccentricity (r p /r e = 0.45). 




Fig. 8. The time-averaged flux radiated by a thin accretion disk 
around rotating neutron and quark stars, with the same central 
density (p c = 10 15 g/cm 3 ) and eccentricity (r p /r e - 0.45). 



4.3. Accretion disks onto neutron and quark stars with fixed 
central density and eccentricity 

As a third example of the comparison between the radiation 
properties of the thin accretion disks onto neutron and quark 
stars we consider the case in which the rotating central compact 
objects have the same central density and eccentricity r p /r e , re- 
spectively, where r p is the polar radius and r e is the equatorial 
radius of the star. The basic fundamental parameters of the rotat- 
ing neutron stars and quark stars are presented in Table [3] 

The physical parameters of the rotating stars have been ob- 
tained by choosing for all the considered equations of state 
a central density p c — 1 x 10 15 erg/cm 3 , and an eccentricity 
r p /r e = 0.45, respectively. We present the radial distribution of 
the energy flux F(r) in this stellar configuration for the neutron 
star models DH, RMF soft/stiff, STOS T = 0.0, T = 0.5 MeV 
and T = 1 .0 MeV, BBBAV14, BBBParis, APR, and for the quark 
star models Q and CFL300, respectively, in Fig. [8] 

Fig. [8] shows that different equations of state with the same 
central density and eccentricity produce rather different geome- 
tries around the compact central object, and in turn determine 
rather distinct emissivity properties of the accretion disks. The 
distribution of the radial flux profiles is similar to the Keplerian 
case, shown in Fig. [5] but with some significant differences. The 



differences between the amplitudes of the disks radiation for the 
central objects with the same p c and r e /r p ratio are smaller as 
compared to the ones obtained in Fig. [5] This is due to the con- 
siderable decrease of the maximal flux values for the models DH, 
BBBVA14 and BBBParis, respectively, as compared to the max- 
ima reached in the Keplerian case. The other EOS types produce 
more or less the same, or somewhat higher, flux maxima as com- 
pared to the case of the Keplerian rotation. By considering the in- 
ner edges of the disk, we find that the radial distribution of their 
location is very similar to the one seen in Fig. [5] Nonetheless, 
the differences between the flux profiles for the STOS models 
are bigger, and there is also a shift of the inner edge of the disk 
to higher radii for the models DH, BBBVA14 and BBBParis, re- 
spectively, which determines a significant reduction in their flux 
values. These features can also be seen in Fig. [9] where we dis- 
play the temperature distribution of the disk surface in thermo- 
dynamical equilibrium. 

Important differences are also present in the disk spectra. 
The black-body spectrum of the equilibrium thermal radiation 
emitted from the accretion disk for each type of star is shown 
in Fig. [10] The quark stars together with the STOS type neutron 
stars produce greater maximal amplitudes of the disk emission 
spectra as compared to the other neutron stars in this group. The 
inverse proportionality between the temperature of the star and 
the maximal amplitudes of the spectra holds again for the STOS 
models, similarly to the Keplerian case. However, the spectrum 
of the quark star with EOS Q and the spectrum of the neutron 
star with the RMF stiff type EOS are decoupled in this case, 
and the disk spectrum for the quark star is more similar to the 
spectrum of the model STOS 1 . The latter has a slightly lower 
cut-off frequency, and these two types of EOS could be discrim- 
inated from each other by comparing the radial distributions of 
the flux and of the temperature of their disk. The spectra of the 
group consisting of the models DH, RMF soft, BBBAV14 and 
BBBParis are similar. In fact, the spectra coincide for the latter 
two models. For this case the different spatial distributions of 
the disk radiation can help again in the discrimination. There is 
no such a problem with the APR type EOS, since it produces 
somewhat higher spectral amplitudes as compared to the other 
cases. The maxima of all the spectra are located in a very narrow 
frequency range between 2 x 10 16 and 2.2 x 10 16 Hz, but the nu- 
merical values of the flux maxima are rather different, more or 
less following those of the flux maxima. 
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Table 3. Physical parameters of the rapidly rotating compact stars with same central density and eccentricity. 
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Fig. 10. The emission spectrum of the accretion disk around ro- 
tating neutron and quark stars, with the same central density 
(p c = 10 15 g/cm 3 ) and eccentricity (r p /r e = 0.45). 

5. Discussions and final remarks 

In the present paper we have considered the basic physical prop- 
erties of matter forming a thin accretion disc around rapidly ro- 
tating neutron and quark stars. The physical parameters of the 
disc - effective potential, flux and emission spectrum profiles - 
have been explicitly obtained for several equations of state of 
the neutron matter and for two types of quark stars, respectively. 
All the astrophysical quantities, related to the observable prop- 
erties of the accretion disc, can be obtained from the metric of 
the central compact object. Due to the differences in the space- 
time structure, the quark stars present some very important dif- 
ferences with respect to the disc properties, as compared to the 
neutron stars. This procedure allows not only to discriminate be- 
tween neutron and quark stars, but also gives some distinct sig- 
natures to differentiate between neutron stars described by dif- 
ferent equations of state, and to distinguish between the different 
phases of the quark matter. 

Another observable physical parameter of the disk that could 
help distinguishing between different classes of compact objects 
is the efficiency e of the conversion of the accreting mass into 
radiation, given by Eqs. (jT7j and Eq. ( fTST l. The values of e mea- 
sure how efficient the energy generating mechanism of the mass 
accretion is. The amount of energy released by the matter leaving 



the marginally stable orbit, or the inner edge of the disk touch- 
ing the surface of the star, and being transferred to the star, is 
the binding energy E ms or E e . We present the efficiency for the 
energy conversion for the different equations of state considered 
in the present paper in Table|4] The first two lines of Table|4]con- 
tain the conversion efficiency of the rotating compact stars with 
the same total mass and angular velocity. 

For this case the values of e are in the range 4% and 8%, in- 
dicating that in this case there is no significant difference in the 
efficiency of the mass-to-radiation conversation mechanism. The 
less efficient conversion mechanism is provided by the STOS 
models, especially at the zero temperature of the star, whereas 
the APR type EOS provides the most efficient conversion mech- 
anism with e = 7.4%. The quark stars with e = 6.8% are also 
efficient engines for the conversion of the accreted mass into out- 
going radiation. 

The values of e for the compact general relativistic objects in 
Keplerian rotation are given in the second two lines of Table [4] 
The case of stars with the same total mass and angular velocity 
is represented in (a), stars in Keplerian rotation are described in 
(b), while the case of the stars with the same central density and 
eccentricity is shown in (c), respectively. 

In this case the quark star model CFL300 provides the most 
efficient mechanism for the conversion, with e = 11.7%. The 
quark stars described by the Q EOS have also a high value of the 
conversion efficiency, with e = 8.15%. The conversion efficiency 
for the models STOS are much higher than the ones in the pre- 
vious stellar configuration. Their efficiencies can reach for this 
case a value of e = 9%. The other models have efficiencies of 
the order of 5-6%. 

The last two lines in Table [4] contains the values of the effi- 
ciency for stellar models with the same eccentricity and central 
density, respectively. The range of the values of e is similar to 
the previous case, although in general the numerical values are 
somewhat lower. We obtain again the highest efficiency for the 
model CFL300, and the smallest ones for the neutron stars DH, 
RMF soft, BBBAV14 and BBBParis, respectively, with e about 
5%. The STOS type neutron stars and the Q type quark stars are 
more efficient, with e = 1 — 8%. 

As shown by the flux integral in Eq. (J6]l, and the explicit 
expressions of the specific energy, specific angular momentum 
and angular velocity given by Eqs. (fTTT i. (fl~2T > and (fT3l , the rather 
different characteristics of the radial flux distribution over the 
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Table 4. The radius of the inner disk edge r,„ and the efficiency e. 



accretion disk, the disk spectra and the conversion efficiency are 
due to the differences between the metric potentials of the neu- 
tron and of the quark stars, respectively. Even if the total mass 
and the angular velocity are the same for each type of the rotating 
central object, producing similar values of £2, E and L, the radia- 
tion properties of the accretion disks around these objects exhibit 
discernible differences. The reason is that the proper volume, and 
in turn the function yf—g, used in the calculation of the flux in- 
tegral, are highly dependent on the behavior of the metric com- 
ponent g n = {drldr) 2 gTf. The latter contains the r-derivatives 
of the metric functions p(r) and y(r), via the coordinate trans- 
formation between 7 and r, which are extremely sensitive to the 
slope of p(r) and y(r). As a result, though the inner edges of the 
disks are located at almost the same radii, the amplitudes of the 
energy fluxes emerging form the disk surface and propagating in 
any solid angle may have considerable differences for different 
equations of state. These features also give the distinctive fea- 
tures in the disk spectra for the various types of central stars. 

In conclusion, the observational study of the thin accretion 
disks around rapidly rotating compact objects can provide a 
powerful tool in distinguishing between standard neutron stars 
and stars with exotic equations of state, that have underwent, for 
example, a phase transition from the neutron phase to a quark 
phase, as well as for discriminating between the different equa- 
tions of state of the dense matter. 
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